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Abstract 

We consider four scalar mesons scattering in large-W QCD. Using 
the worldline formalism we show that the scattering amplitude can be 
written as a formal sum over Wilson loops. The AdS/CFT correspon¬ 
dence maps this sum into a sum over string worldsheets in a confining 
background. We then argue that for well separated mesons the snm 
is dominated by flat space configurations. Under additional assump¬ 
tions about the dual string path integral we obtain the Veneziano 
amplitude. 


The 1968 seminal paper by Veneziano marks the birth of string theory 
(“dual resonance model”). It concerns a proposal for the scattering amplitude 
of four scalar mesons, as follows [1] 


A{ki, k2, ks, ki) 


r(— q;(s) — a{t)) 


( 1 ) 


with q;(s) = q;(0) + a's, a{t) = q;(0) + a't and s, t are Mandelstam variables. 

While the amplitude (1) is appealing due to its theoretical and phe¬ 
nomenological properties, it lacks a derivation from Quantum Chromody¬ 
namics (QCD). It is important to clarify the exact relation between the 
above amplitude and QCD. In particular, we wish to hnd whether there ex¬ 
ists an approximation where we can derive the amplitude (1) from QCD. The 
topology of the string diagram (“disk amplitude”) suggests that it should be 
related to large-W QCD with hxed Nf, namely to the’t Hooft limit of QCD 
where quark loops are suppressed (quarks are “quenched”). As we shall see, 
the route from QCD to (1) involves additional approximations which are not 
controlled by QCD parameters. In particular, we assume that Wilson loops 
are calculated by flat space string worldsheets. 

The hrst attempt to relate the Veneziano amplitude to held theory was 
made in 1970 [2, 3], even before the birth of QCD. Sakita and Virasoro argued 
[3], using a scalar held theory as a model of strong interactions, that the held 
theory amplitude should look like a dense hshnet, see hg.(l). Moreover, they 
argued that at high orders in perturbation theory the gaps between the holes 
in the hshnet close and the amplitude should resemble a string amplitude, 
namely the value of the amplitude should correspond to the area of the string 
worldsheet. While this idea can easily be implemented in large-W QCD (with 
a hshnet made of gluons), it cannot be the full story as it lacks the essential 
non-perturbative ehects which are responsible for conhnement with a string 
tension and a mass gap. It is therefore necessary to address the problem by 
using a formalism that incorporates non-perturbative ehects. 

We consider SU{Nc) QCD with Nf massless quarks in the’t Hooft limit, 
where g'^N^ and Nf are kept hxed. The partition function takes the form 


Z = j DA^expi-Syu) (det(* P)ff 


( 2 ) 


Let us write the fermionic determinant as a sum over (super-)Wilson loops 
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Figure 1: Perturbative description of the four mesons scattering amplitude as a “fishnet” 
diagram in planar QCD. 


by using the worldline formalism [4] 


with 
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1 dT 

2)0 ^ 


X 


J Vx^Vtl)^ exp 


exp(iV/r[A^]), 



—x^x^ + 



(3) 


X TrPexp |z ^ dr (^A^x^ - } , (4) 

where p, i/ = 0,3. 

The expansion of the exponent in powers oi Nf 

exp(iV/r) ~ 1 + NfV + 0{N‘j) (5) 

is in fact an expansion in Nf/Nc, since (hFihF 2 ---Wd)conn. N^NlIKf. 
Therefore, as we shall see, in the’t Hooft limit the dominant class of diagrams 
for meson scattering involves one Wilson loop [5]. 

Let us consider the scattering amplitude of four scalar mesons 


A{xi, X2, Xs, X4) = {qq{xi)qq{x2)qq{x3)qq{x4)) . ( 6 ) 


For simplicity we assume flavor-singlet mesons. Flavored mesons can be 
easily incorporated by adding “Chan-Paton” factors, resulting in an overall 
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group theoretic factor of Tr (A^A^A^A'^). In order to calculate the amplitude 
we add to the action a source of the form f d^x J{x)qq 



(7) 


and differentiate the partition function with respect to J{xi), J{x 2 ), J^x^) 


6 5 5 5 


{qqixi)qq{x2)qq{x3)qq{x4))) 


\ogZj\j=o (8) 


6 J{xi) 5 J{x 2 ) 6 J{x 3 ) 5J{xi) 


In the worldline formalism, such a variation excludes Wilson loops whose 
contours do not pass through {xi, 0 : 2 , 0 : 3 ,X 4 } [ 6 ][5]. 

Thus in the’t Hooft limit the leading large-W expression for the scatter¬ 
ing amplitude can schematically be written as 



, X2, X3, X 4 ))yM • 

(9) 


In order to simplify the notation we omitted the worldline fermions in the 
above expression. The Wilson loops are calculated in the large-iV pure Yang- 
Mills theory. The sum is over all sizes and shapes of Wilson loops that pass 
through the points {xi, X2, X3, X4}. The QCD amplitude admits the topology 
of a disk and it resembles the string disk amplitude (2) [7]. 



Figure 2: Worldline description of the four mesons scattering amplitude. It resembles 
the string disk amplitude. 

In order to proceed we use holography to evaluate (fF(xi, X 2 , X 3 , X 4 ))ym 
( instead we could have used the lattice strong coupling expansion). The 


3 





holographic prescription for evaluating the Wilson loop expectation value is 
to hud the minimal worldsheet that ends on the boundary of the space (anti 
de Sitter in the case of A/” = 4 super Yang-Mills theory) and passes through 
{Xi,X2,X3,X4} [ 8 , 9 ]. 

We therefore propose that the field theory path integral over all sizes 
and shapes of contours that pass through {xi, 0 : 2 , xs, 0 : 4 }, translates into the 
Polyakov path integral for string worldsheets with a disk topology that end 
on the boundary of the space and pass through {xi, 0 : 2 , 0 : 3 , 0 : 4 } 

^ J Vx^^ir) exp (^- J dr^xf^ {W{xi,X 2 ,x^,xa)) = ( 10 ) 

j Vg°‘^Vx^ eyrp {^ 1 ^, j d^a y/gg^'^d^x^dfiX^GMN^ \{xi,x2,x3,xi} 

with M, N = 0,D [D = 9 for the superstring). The above proposal (10) is 
conjectured to hold for all gauge/gravity dual pairs. The information about 
the specihc gauge theory is encoded in the metric Gmn- We will use Witten’s 
model of back-reacted W D4 branes compactihed on a thermal circle [10]. 
The conjecture is [10] that type IIA string theory on 

rirG 

= {UIR)^l^{r]^,dxi‘dx‘‘ + f(U)dT^) + + U^dttl) 

f(U) = 1 - , (11) 

with anti-periodic condition for NS-R modes on r is dual to large-Y Yang- 
Mills theory. Unfortunately, in this model the scale of the Kaluza-Klein 
modes and the scale of conhnement are of the same order, S 
[11]. Note that since we use the type IIA superstring, we should use in 
(10) the Polyakov action for the superstring, namely we should incorporate 
worldsheet fermions. 

A calculation of Wilson loop expectation value in Witten’s background 
yields an area law for large loops, (W) = exp(—SA), with a string tension 
^ reason is that for a sufficiently large loop the string 

will quickly drop to U = Ukk and the main contribution to the action is 
essentially by a flat space worldsheet that resides close to Ukk [ 12 ]. 

It is therefore anticipated that if the points {xi, 0 : 2 , 0 : 3 , 0 : 4 } are well sepa¬ 
rated from each other \xi — Xj\ 3> l/\/S, the sum over string worldsheets will 
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be dominated by flat space conflgurations, as depicted in flg.(3). This claim 
is supported by the analysis of ref. [6], where the authors showed that under 
the assumption of confinement (area law) the sum over Wilson loops is domi¬ 
nated by a saddle-point that corresponds to the IR piece of the worldsheet in 
fig. (3). Our claim is indeed similar to that of ref. [6]: assuming that Wilson 
loop expectation values are calculated by a flat space string configuration 
implies that they admit an area law. 


U 


IR cut-off 



Figure 3: A string worldsheet in a confining background. The string “sits” at the IR 
cut-off. 

String worldsheets that span in the compact directions (r and the four- 
sphere) are expected to admit a larger action than string worldsheets that 
reside in the flat four dimensional space. 

If this is indeed the case, we proceed as follows. We first fix the gauge 
by choosing a flat worldsheet metric and add ghosts. At tree level the ghosts 
decouple from and therefore will not play a role in the calculation of 
the disk amplitude, so they will be omitted. Since we anticipate that the 
contribution of compact directions (r, and their superpartners) is not 
dominant, we assume that we can carry out the integration over these fields as 
if the metric is flat. In fact, if the metric was U independent, the integration 
over the compact directions becomes indeed trivial. 

Under these assumptions the calculation is crudely approximated by the 
path integral 

j Vx^^VU exp j (i^(T{{d^x^^d°‘x''G^y) + {dJJd°‘UGuu)}^ I{xi,x2,x3,x4}, 

( 12 ) 
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with the 5d metric 


ds^ = {U/Rf^\,dx^dx^ + {R/Uf^{j^). (13) 

Due to the metric singularity, the path integral (12) is dominated by U in 
the vicinity of Ukk- This is the case for large worldsheets, where the string 
“sits” at the horizon. Thus, let us insert 5{U — {Ukk + e)) into the path 
integral, suppressing quantum fluctuations in the U directions. Later on we 
will discuss what is expected to happen if we omit the delta fnnction and 
allow flnctuations in the U direction. The modihed path integral reads 

j Vx^^VU5{U - Ukk - e) x 


X exp j d^a{{daX^d'"x''Gf,^) + {de.Ud'^UGuu)}^ |{a;i,a; 2 ,a; 3 ,:r 4 } 



with G^u the flat 4d metric 

ds"^ = (Ukk/R)^^‘^V tJ.ydx^dx’^ . (15) 

It is easier to calculate the amplitude in momenta space. Let us carry out 
the computation with external states of kf = 0. This can be justihed by as¬ 
suming a scattering of fonr massless scalars. A physical case that corresponds 
to massless (pseudo-)scalars is the scattering of pions - the Nambu-Goldstone 
bosons that correspond to chiral symmetry breaking. 

To this end we assnme that the appropriate scalar vertex operator to 
insert into the path integral takes the form 

V j dy exp{ik ■ x) (16) 

(the vertex operators sit on the disk boundary). Note that due to our as¬ 
sumption that U = Ukk we can use flat space vertex operators. Since our 
bulk theory is type IIA the vertex operators involved in scalar scattering are 
anticipated to be snpersymmetric and not the assertion (16). The proposal 
(16) is motivated by the fact that snpersymmetry is broken. Since the only 
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massless states in the gauge theory are the massless Nambu-Goldstone bosons 
(the “pions”), the vertex operator (16) of the dual string must correspond 
to these states, hence it is suggested to be a purely bosonic vertex operator. 
The assertion (16) is the simplest operator which involves massless bosonic 
helds. 

Inserting vertex operators in the path integral is not equivalent to a simple 
Fourier transform, but it rather leads to an element of the scattering matrix, 
see [13] for a recent discussion. 


A{ki, k2, ks, ki) = 


(17) 


fvx^ n /%.exp(.Ml/*))exp(-S 


d‘^a daX^d°‘Xfj, . , 


with yi a coordinate on the disk boundary. The scattering amplitude is 

_ \2a'{RIUKKf/^kj-k, 


A{ki, k2, k3, ki) ~ (5(^ ki) / JJ dyi n(» -yi) 

i i j<l 


(18) 


with i,j,l = 1..4. Since we started with type IIA string theory and we 
carried out a calculation in flat space, the result is the celebrated Veneziano 
amplitude for the superstring with a string tension S. 

An important issue is the suppression of fluctuations in the U direction. 
Although classically the string sits aX U = Ukk, quantum mechanically 
it fluctuates. Incorporating such fluctuations is not an easy task, since it 
involves a solution of the non-linear sigma model. The fluctuations near 
Ukk are expected to affect the low lying spectrum of the string and to lead to 
deviations from the exact linearity of the Regge trajectories that characterise 
the Veneziano amplitude. 

In conclusion, we have shown that the large-A^^c, hxed Nf, QCD disk 
amplitude that represents scalar mesons scattering can be written as a sum 
over Wilson loops. By using holography the calculation is mapped into a sum 
over string worldsheets. Carrying out this sum for well separated mesons (or 
small Aki, corresponding to the low-energy regime) ignoring fluctuations in 
the U direction, leads to a flat space calculation which yields the Veneziano 
amplitude for the type IIA superstring. 

Note that while we used holography to relate Wilson loops in the gauge 
theory side to string worldsheets, we did not need to use the intuitive rela¬ 
tion between mesons and open strings. Our approach is different from other 
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approaches that relate gauge theory scattering amplitudes with string am¬ 
plitudes. In particular, in contrast to Polchinski and Strassler [14], we are 
interested in meson scattering and not in glueball (pomeron) scattering and, 
moreover, we use the method of [5] in order to introduce mesons as opposed 
to the more common method of using a probe brane. 

We therefore demonstrated how string disk amplitudes are related to 
QCD amplitudes. Our analysis can be easily generalised to the case of k- 
mesons scattering. A more challenging task is to take into account the space- 
time curvature and fluctnations in the U directions in order to obtain cor¬ 
rections to the Veneziano amplitude. If there were a gravity dual to large- 
QCD and if we could carry out the full path integral we should have obtained 
an amplitude that describes correctly the DIS regime of QCD. Finally, it is 
expected that some of the perimeter effects will account for corrections due 
to quark masses, as in [15]. 
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